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Abstract 



^><-^ I The framework for soft leptogenesis minimally extended with a DM sector is studied. 

A heavy singlet neutrino superfield acts as the source for (s)lepton asymmetry and by 
coupling to the singlet DM superfield it produces a DM particle density through decays. 



pi I The nature of DM generated is twofold depending on whether the Yukawa and DM 

couplings are either small or large. With sufficiently small Yukawa and DM couplings DM 
annihilations into MSSM particles are slow and as a consequence all DM particles form 
the DM component. The solutions to Boltzmann equations are given and the dependence 



>: 

■ between the DM masses and coupling are presented in this weak coupling regime. Also, 

m ■ 

ff^ ■ the behavior of the efficiency of producing asymmetric DM is determined with weak 



couplings. We note that a different outcome arises if the couplings are larger because 
then the ADM component is dominant due to the effectiveness of DM decays into the 
MSSM sector. 



u ' 1 Introduction 



Our current understanding of the universe is that around 5 % of the energy density belongs to 
ordinary baryonic matter while the rest stems from dark matter (DM) and dark energy with 
proportions 22 % and 73 %, respectively [Ij. The existence of all three components is more 
or less a mystery and points towards physics beyond the Standard Model (SM). The fact that 
there seems to be almost exclusively matter and no antimatter in the universe suggests there 
is matter-antimatter asymmetry. Since the SM cannot accommodate sufficient baryogenesis 
on its own, baryogenesis mechanisms have been constructed in various Grand Unified Theory 
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(GUT) scenarios. Alongside GUT models, baryogenesis via leptogenesis [2] has been studied 
extensively over the last twenty years or so. 

As for DM, a viable candidate arises in the Minimal Supersymmetric Standard Model 
(MSSM), namely the Lightest Supersymmetric Particle (LSP) [3l H]. Usually, the lightest 
neutralino (i.e. a combination of the superpartners of the neutral gauge bosons and higgs 
sparticles), or in supergravity theories the gravitino, assumes the role of the LSP, and the 
success of the LSP in producing the correct DM abundance depends on the adjustment of 
parameters within the MSSM. However, since the DM and baryonic energy densities are close 
to each other, i^Divi/^b ~ 5, it is natural to assume that DM and baryogenesis are somehow 
related. The connection between baryogenesis and DM has also been widely explored [S]. 
Moreover, the possibility of the DM energy density being the product of an asymmetry has 
led to the concept of asymmetric DM (ADM) where DM and DM antiparticles are distinct 
from another and produce an asymmetry analogous to baryons. Note, however, that as ADM 
scenarios incorporate MSSM, also LSP is always included in the ADM models but is not 
necessarily good candidate for DM. 

Over the years various realizations of leptogenesis have been devised [6]- [28]. They are 
usually motivated by SO (10) scenarios where B — L is a gauge symmetry, and it is natural to 
have a heavy right-handed SM singlet neutrino that can also be responsible for light neutrino 
masses via the seesaw mechanism [2H 125]. Adding a heavy singlet neutrino superfield into the 
Minimal Supersymmetric SM (MSSM) gives rise to soft leptogenesis where the soft breaking of 
supersymmetry allows the scalar superpartner of the heavy neutrino to generate the (s)lepton 
asymmetry f2E[ W \. In contrast to non-SUSY models, only one family of heavy singlets is 
required. 

Furthermore, various scenarios incorporating leptogenesis and DM have been introduced 
[29]-j37]. In this article, we investigate a model that encompasses soft leptogenesis and DM 
generation from the same heavy singlet superfield. Previously a model with SU{2)l triplet 
DM and supersymmetric leptogenesis has been studied [2B]- Our model is a minimal extension 
to the Minimal Supersymmetric Standard Model (MSSM) and it consists of MSSM augmented 
with a single SM singlet DM superfield and a heavy singlet neutrino superfield that enables 
soft leptogenesis after soft SUSY breaking. In section |5] we introduce our model and present 
the CP violation parameters, in section [2] we solve the Boltzmann equations and the resulting 
baryon and DM abundances, and in section 2] we present our conclusions. 
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2 The model 



Models incorporating leptogenesis and DM generation have the potential to produce the ob- 
served matter and DM energy densities and offer a possible explanation as to why the abun- 
dances are similar in magnitude. So far the MSSM has been considered as one of the most 
appealing extension to the SM, and with the addition of a heavy singlet (s) neutrino, soft lep- 
togenesis can arise. To account for DM, the soft leptogenesis framework can be extended with 
an additional DM sector. We have augmented the MSSM with a singlet superfield consisting of 
a heavy neutrino and sneutrino and a singlet DM superfield. The heavy (s)neutrino couples to 
the MSSM (s)leptons, Higgses and higgsinos, which gives rise to the soft leptogenesis scenario. 
The DM sector couples to the heavy (s)neutrino as well and this provides additional decay 
channels for the singlet (s)neutrino and a source for DM production through these decays. The 
superpotential is 



(2.1^ 



The discrete symmetry Z2 imposed on $ forbids terms ~ $^ and ~ NN^ . After adding 
soft supersymmetry breaking terms we get the part of the Lagrangian that includes the soft 
leptogenesis and darkgenesis interactions 



C = X(l)Nx + VaHuL^N + yaL^HuN + h.c. 

+ ^y^M*e''^/'L^H^ + h,e-''^/'<P^ + h^e-^'^/'l^H^ + h.c. 



-^y^M*e'^''L^H^ + -^cie-^^/V' + ^c^e^^^/^L^if. + h.c. 
where the sneutrino mass eigenstates are 

V2i V 



(2.2) 

(2.3) 
(2.4) 



The DM particles x cind represent the DM fermion and scalar, respectively. 

The lepton and DM asymmetries arise from the self-energy loop diagrams shown in Figs. 



1(a) 1(f) While the process A''-!- — )■ xx'^ does not produce a DM asymmetry, i.e. an asymmetry 
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between x and x*^, the decays N± — >■ 00 and N± — >■ 0*0* can create an asymmetry between 
and 0*. This is the source for possible ADM in our model and the asymmetric component 
consists of scalar DM. The CP violation parameter for (s)lepton production is defined as 



[r(iv+ ^ /) - r(7v+ ^ /) + r(7v_ ^ /) - r(7v_ ^ /)j 
E/ [r(iv+ ^ /) + r(7v+ ^ /) + r(7V_ ^ /) + r(7V_ ^ /)] ' 



(2.5) 



where / stands for the fermionic {laH^) and bosonic {laH^) final states. The DM asymmetry 
is defined in a similar manner: 



I (^N+ ^ 00^ - r (^N+ 0*0*) + Ti^N_^ 00) - t(^N_^ 



r (^A^+ ^ 00) + r (^N+ 0*0*) + r (^A^_ ^ 00) + r (^a^_ ^ 0*0*) 



(2.6) 
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The temperature-dependent CP violation parameter for the (s)leptons is 

4 |M| 



167r{Ml^-Miy + Ul_ vl 
|M|2 \A\' 



2\A\\M\ cos (eA + ^ + eM- Ob) ^meff 



+ |A| |M| |ci| COS (-^A + I + - ^6 + ^ci) ) 



IGttM 



M% -Ml 



|M|2 |M| 



167r(Ml_-Mj^)2 + n^, 
X (^2|M||A| cos (Om + OA-Ob + ^) 

+ |ci||A||M|cOS (^ci-^A + ^M-^ft 



-meff 



2/7 167rMjy_ 



Ml -Ml 



\M\m,s f 



|M| 



167r(Ml_-Mj^)2 + n^^ ^2 2 V 

+ |ci||A||M|sin(^A-^c,-^^M + ^^fe 
4 iMlmef, 



2|M| I A| Sin (-^M - + ^ft) V^^eff 



167r ^- 



167r(M2 -M2 )2 + ni_ ^2 2 



|M| 



2| A| |M| Sin {Ob -9a- ^m) Vmeg 



+ |ci||M||A| sin - e,, +6^- Om) J^iv+J^ 



Cf 



X 



(|M|2 + + 2\M\ \A\ cos (^M -e,-r 9 a)) ^^eff^ 



+ (|M|2 + - 2|M| 1^1 cos (^M - ^6 + Oa)) ^meff '''' 



87rM^_ 



1^1 ,^ |M| cj. 



(2.7) 



where we see the influence of the DM sector couphngs in the N± self-energies and N± — > 
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transitions. We have used C2 = A Va and rUeS = v'^ Yla IVal'^/^- The net asymmetry is 

4 '^'~N^ -'^'~N /|M|2|A|2 



Ml -Ml 
IQTiiMi -Ml )2 + n2 V 8 



Cl 



2 / IQirMfi 



X 2|v4||M|cos 



— + — t/b I — 



+\c,\\M\Xsm{et-e,, + ex-eM 

4 ^l_-^L /|ci|2 lAplMp 



16n{Ml -M| )2 + n 



167rMjv_ 



(^2|M||A| cos (^M + - ^fe + ^) ^meff 
+ |ci||A||M|sin(eA + ^^;,-^c, -^^A/) 



X 



C0 



with the factors [T2 



CP = (l-x,-x^)^(l + x,-X5)'-4x,(l-/n 
= 7(l + XH-a;^)'-4a;^,(l + /l,^) (l + ) , 

C0 = v/r^(i+/r)(i+/r)' 



M2 

where the thermal masses are 

mi{Tf 



T 
1 



. 9y 3^ 
92 + 4 ' ' 



IT' ( 3 2 9y 3 
2 V8^^ 8 4 ' ' 



2/39 



9l 



T 

2 \8^^ 8 



ml + m|,„ft + A^Tl 



(2.8) 



(2.9) 

(2.10) 
(2.11) 

(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 
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The phases 6m and 9b can be rotated away and are set to zero and so M is real and we can 
replace \M\ — )■ M. The factors cf,b,<i> reflect the production thresholds of MSSM particles 
and DM boson (p imposed by energy conservation. Since we normalize the DM asymmetry 
to the DM boson production, the factors cancel and consequently there is no temperature 
dependence in ^dm- We have also kept the soft terms at the vertices for completeness, usually 
these are neglected [27j . 
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3 Boltzmann equations 



We move on to solve the relic baryon and DM abundances this scenario generates. We assume 
the Yukawa couplings rries and DM coupling |A| are sufficiently small and consider only pro- 
cesses with cross sections 0{y^) and 0{\^). Due to this, DM annihilation processes are left 
out, which rules out ADM as the primary source of the DM energy density. Thus, we write 
the Boltzmann equations also for the total abundances of x and 0. The particle abundances 
are expressed in terms of the number densities normalized to entropy density Yi = rii/s and 
their evolution with changing z = M/T is to be determined. Our Boltzmann equations are 

dYjv 



dz 



N 



dz 



dz 



dY, 



dz 



dz 



N 



yeqx 
^ N I 



(3) 



AS, 



N 



+ Df + 



e,{T) (y^ - 2Yp 



D 



N^eq 



Y^q I 



2 



SI 



(8)' 



" N 



■ N N 



(4) 
t 



+Yn (2Sf^ + ^) + 2F^^ (Sf ) + 5f ) + si 
+2Y1^ (^f ) + SP + SP) + S,2 (^2Yp + ^ 



D 



Y, 



yeq ^ 9 -^JV - ' 



■ N 



Ami 

-J^Dnxx+\I^ + 



{ml 



mj 



"ml + ml 



M4 



M2 



-D 



Nx4' 



(3.1) 



(3.2) 



(3.3) 
(3.4) 



dY^ 
dz 



X - 



Am\ / ( 



mt 



M4 



M2 



Nx<p 



(3.5) 



X - ^yp 



(3.6) 



with Yp ^ Yp^ and Y^/2 ^ [39j. Recall that we take into account 0{y^) and 0{X^) 
decay rates and scattering diagrams involving the (s)top. The expressions for the decay D 
and (s)top scattering rates S can be found in appendix A. The solutions to these Boltzmann 
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equations are shown in Figs. [2llU The baryon abundance is related to the (s)lepton abundance 
through 

Yb = -^>AL.... (3.7) 

We have chosen m^Q = 10'^* TeV, M = 10^ TeV, Mh=125 GeV, \A\ = 1 TeV, |A| = 10"®, 
|ci| = 0.001 TeV and m$ = 0.1 TeV0 The evolution of the Yb abundance is similar to those 
found in [271 [39]. 

The DM abundances and are not sensitive to changes in the DM mass parameter m$ 
as long as it is a lot smaller than the heavy scale M associated with the singlet (s)neutrino, 
m<s, -C M. It is the coupling A and the positive mass dimension soft coupling ci that largely 
determine the order of magnitude in the DM abundances. With suitable choices of these 
couplings, the abundances Y^ and Y^ attain final values ~ 10^^^ suggesting that both the 
scalar and fermionic DM particles can be viable DM candidates. 

At small couplings A ~ 10~^ and ci ~ 0.001 TeV, the asymmetric component is extremely 
tiny with abundance ~ 10^^'' and does not play a role in DM energy density. However, an 
increase in both DM sector couplings causes the magnitude of to rise and also compels to 
consider scatterings where DM particles are destroyed and MSSM particles are produced. So, 
in order to have ADM as the primary source of all DM, it is crucial to find a balance between 
sufficient DM production through asymmetric decays and scattering processes 

XX^ LaHu, XX^ LaH^, 
X(p L^H^, x<P L^Hl, 

where x'^ denotes the right-handed fermion obtained from x through x" = (x^s)*? where 
_ (^_^^2^"/3^ believe it is possible to create conditions where the DM particles annihilate 
and the asymmetric component remains and forms a part of final DM abundance although 
this scenario entails a different parameter space in terms of DM couplings and other parameter 
inputs. We are studying this alternative in an upcoming paper. 

We have explored the allowed parameter space and made scatter plots of the DM masses 
versus the DM coupling |A|, Figs. [5ll6l The fixed parameters are meg = 10~^^ TeV, M = 10^ 

^If the heavy (s)neutrinos are supposed to be produced thermaUy, large M imphes reheating temperature 
Tji > M . In many inflatory supergravity models this high Tb. leads to overproduction of gravitinos. However, 
recently it has been proposed that with the present value of the Higgs mass about 125 GeV [42J the reheating 
temperature constraint could be relaxed, Tr ^ 10^ — lO^*' GeV, not in conflict with our scenario jiB] , 
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Table 1: The DM sector parameters and their scanned values. 



Parameters 


Scanned values 


\M 


10-10 _ 10-8 


|ci| 


10-5 - 10-3 TeV 




0.01-1 TeV 


2 

"^*soft 


10-^-1 TeV2 



TeV and A = 1 TeV while |A|, |ci|, m<i, and m$soft are varied as shown in Table 1. The 
range for the soft parameter |ci| has been chosen so that |ci| ~ |A|M. The DM coupling 
10^^° < |A| < 10^8 is the most favorable region of producing the correct DM abundance which 
is why this particular region is studied. The DM mass parameters are chosen to be in the 
weak scale and even in the 10 GeV range which brings the typical ADM regime closer. The 
required DM abundances come from the relation 



K\+m,\Y,\<l^T (3.8) 



and the baryon abundance has to satisfy 

IYbI <8.8x 10-" (3.9) 

which we have used as constraints in extracting the allowed parameter space. Known parameter 
values from observations are T = 2.725 K, Q^m = 0.11//i^ and = 2.47 x 10^^ /h"^. We have 
taken the lower limit in (13.91) to be 30 % of the observed baryon abundance while in (13. 8p the 
lower limit is 90 % of the observed value. Figs. [SHH] reflect the hyperbola-like boundary set by 
Eq. (13. 8p . and there is an indication of slight preference of |A| values in the ~ 10^^ ballpark. 
A similar observation of the favorable magnitude of |A| can be made from a plot with the 
left-hand side of (13. 8 p plotted against |A|. 

We have also determined the dependence of efficiency factors related to baryon and ADM 
production on the Yukawa and DM couplings. The efficiency factors t]b and r^oM are shown in 
Figs. and they are given by 



Vb,dm 



A<t> 



2£l,dm1^~^(^ ~ 0) 



(3.10) 



The ADM production efficiency falls with rising while with increasing |A| //dm increases 
and after reaching a maximum a steep decline ensues. The fall in t^dm is due to the fact that 
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while |A| increases, DM production climbs up but sooner or later washout processes overcome 
the production. This observation is also made in [3S]. Of course a more accurate prediction 
of the DM efficiency factor can be made when the scenario better fulfills the criteria for the 
existence of ADM, which would mean larger DM couplings |A| (and |ci|) and the inclusion of 
DM annihilations into MSSM particles. 




(e) (f) 
Figure 1: The relevant loop diagrams contributing to (s)lepton and DM asymmetries. 

4 Conclusions 

We have studied soft leptogenesis with an additional singlet DM superfield that couples to the 
singlet neutrino superfield. This framework allows the production of DM through the heavy 
(s)neutrino decays. In the weak coupling regime we have neglected suppressed higher order 
scattering processes, which does not make the ADM scenario viable due to the lack of DM 
annihilation channels. The absence of annihilations makes the DM produced in the tree level 
decays N± — )■ xx'^t ^± ~^ 00 cind — > x0 stable. 

It turns out that sufficient levels of baryons and DM can be produced in accordance with 
observations without the need for ADM. This is supported by the evolution of the final abun- 
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0.01 




Figure 2: Baryon \Yb\ and sneutrino abundances with M = 10^ TeV, rries = 10 TeV, 
A=l TeV, A = 10-^ and m$ = 0.1 TeV. 

dances Yg and Y^^^^ and parameter space related to the DM sector. We have produced scatter 
plots depicting the overall dependence of DM masses m^^fj, versus the DM coupling A with 
cosmological observations of the baryonic matter and DM energy densities as constraints, Eqs. 
(13. 8p and (13. 9p . The DM coupling A and its soft counterpart Ci are independent from the 
MSSM and soft leptogenesis parameters and thus their magnitudes are fixed by the present 
DM abundance, and also the baryon abundance through loop effects in El in Eq. (12. 7p . The 
plots in Figs. |5] and E] exhibit the fact that the observed DM energy density restricts the 
allowed region by making a hyperbola-like boundary. 

For future studies, it would be interesting to see what is the correct parameter space that 
allows for ADM in this scenario. We believe this happens once we take into account the 
annihilation processes for the DM particles. The leading annihilation channel for scalar DM 
00 — )■ LaHu comes from the F-term potential. Other channels are mediated by and N±. 
Thus, a stronger coupling of the DM sector can not only lead to higher levels of the asymmetric 
abundance but also makes the DM annihilations significant which is a necessary precursor 
for the existence of ADM. 

Acknowledgments The work of H.K. was supported by Alfred Kordelin Foundation. 
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|i'^l/10'^(-6) 



Figure 3: Scalar DM and sneutrino abundances with M = 10^ TeV, rUes = 10 ^'^ TeV, 
A=l TeV, |A| = 10-^ m$ = 0.1 TeV and |ci| = 0.001 TeV. 

A Decay and scattering rates 

The (s)top scattering rates are given by [101 EHl SI] 



where 



Si'' 



(2) 



Ks 
6 



(4) 



5. 



22 



Ks 
6 

6 

6 

3at 



(7) 



6 



K, 



M w„ 
18af— meff— 



M MmesKi{z) 
87rm{z =1)^ z K2{z) 



K. 



ss 



M16%%5/2y^^;2' 



M2 M1627r%5/2 /^y2 



(A.l) 



(A.2) 
(A.3) 



13 



0.01 




\eq 



|F.|/10'^(-6) 



Figure 4: Fermionic DM \Y^\ and sneutrino abundances with M = 10^ TeV, rries = 10 
TeV, A=l TeV, A = 10-^ m$ = 0.1 TeV and |ci| = 0.001 TeV. 



The reaction rates are found by integrating the functions (x = ip/ z^) 
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Figure 5: Scatter plot of vs |A| with the present baryon and DM abundances as constraints, 
Eqs. (1231) and (1^ . 



where ah = (Mh/M)'^ with Mh standing for the mass of the Higgs boson. The integrals 
determining the reaction rates are of the form 

K2iz) 



The (s) neutrino decay rates are given by 
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Figure 6: Scatter plot of vs |A| with the present baryon and DM abundances as constraints, 
Eqs. (1231) and 
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